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^ Abstract 

f^ We demonstrate the conventional Jarzynski relation (JR) is violated for a non-Markovian process 

^_^ with colored noise. As an example an exactly soluble model is considered with a simple protocol 

O for the external work performed on the system along a non-equilibrium trajectory. For that model 

H we derive an exact expression for the dissipative energy in terms of an arbitrary correlator of the 

Q^ noise characterized by an autocorrelation time tc- As the result we find corrections to the JR in 

terms of tc- In the limiting case of a Gaussian process as well as an infinitely slow process the 

conventional JR is retained. The result is valid for an arbitrary colored noise. 
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I. INTRODUCTION 

During the last decade, initiated by the seemingly remarkable relation derived in pQ, a 
number of exact relations have been derived for non-equilibrium processes. The prominent 
equation in that context is the Jarzynski relation (JR) [1]. This relation joins the change 
of the equilibrium free energy AF = Fj — Fi between two configurations of a system, an 
initial (z) and a final state (/), with the non-equilibrium work function W spent in driving 
the system from the i to the / state: 

(e-/^^) = e-^^^ (1) 

Here (3 = {kBT)~^ is the inverse temperature of the heat bath to that the system is coupled, 
ks is the Boltzmann constant. The relation is a highly non-trivial extension of the ther- 
modynamic inequality W > AF valid for fixed temperature [2]. The substantial difference 
consists of the inclusion of an underlying dynamical process indicated by the average pro- 
cedure in Eq. ([I]). Originally the JR was derived using a deterministic dynamics based on a 
Hamiltonian formulation [1] . Later the relation had been extended to a Langevin dynamics 
[3] or a discrete Markov equation [H |3]. The validity of the JR can be demonstrated for 
an exactly soluble models P|. Recently the JR is also applied in single-molecule pulling 
experiments [71 El E] in order to measure the folding free energies, and in calculations within 
the kinetic theory of gases [TOl [TT]. Likewise, the non-equilibrium dynamics of gene ex- 
pression and its relation to the JR had been studied in [12j. An extension of the proof of 
the JR, solely based on specific equations of motion without any further assumption and 
non-Hamiltonian dynamics, has been proposed in [13]. Further, the relation between super- 
symmetry and non-equilibrium work relations was discussed in [H] . Different aspect of the 
JR had been explored in several other papers such as in [15j, where the entropy production 
along a stochastic trajectory had been analyzed and a generalized fluctuation theorem could 
be derived. Time symmetric fluctuations and their relation to the JR are obtained in [T6] . 
JR in quantum systems had been studied in [T71 [TS] . A special problem discussed in that 
context is to find out an expression for the characteristic function of work performed on a 
quantum system [19]. Work fluctuations in quantum spin chains were studied recently [20] . 
The broad spectrum of papers devoted to the JR and related problems suggest that the JR 
is generally valid. In the present paper we show for a simple example that the conventional 
JR is violated for a non-Markovian process. Especially, we calculate explicitly corrections to 



the JR in terms of a parameter which characterizes the non-Mar kovian behavior. Obviously 
one is confronted with different interpretations concerning the fact what a non-Markovian 
process means. Whereas a history dependent process with a white noise is classified as a 
Markovian process [2lj, sometime a process with memory and white noise is declared as a 
non-Markovian process [22] • In this paper we follow the classification scheme given by van 
Kampen [23]. In that context a process subjected to a colored noise is denoted as a non- 
Markovian process. For such kind of processes we demonstrate the breakdown of the JR. 
As shown directly the reason for the violation of conventional JR consists of the coupling of 
different time scales. In our approach there appears a term, by which the work protocol and 
the noise term are coupled. Whereas in case of a white noise the coupling between noise and 
working function is instantaneously, in case of a colored noise, the protocol is coupled to all 
former time scales of the noise. Eventually this coupling leads to a modification of the JR. 

II. MODEL 

To demonstrate the violation of the Jarzynski relation we consider the simple model of a 
damped oscillator which follows the equation for the variable x{t) 

mx + 70: + mujlx = m[K{t) + r](t)] . (2) 

As usual m is the mass, 7 is a damping parameter and ujq is the frequency of the undamped 
oscillator. The function K{t) is an external force. For simplicity we assume the following 
protocol for that external force 

-00 < t < 

t < t < r (3) 



K{t;T) 



Ko 



T 

Kq t < t < 00 



The last relation means that the system is initially in equilibrium. At time t = the system 
is perturbed by a time dependent force. After a time interval r > the external force is 
switched of and the system tends again to an equilibrium state. Notice that the limiting case 
r -^ 00 characterizes an infinite slow process. During such a quasistatic process the systems 
remains in equilibrium. Aside from an external force the system should be subjected to a 
stochastic force ri{t) with the properties 

(r/(t)) = 0, {v{t)v{t')) = S{t-t'). (4) 



The probability measure corresponding to Eq. Q is given by 

P{r]) = ATexp (- f dt f dt'r]{t)S{t-t')r]{t')\ , (5) 

where A/" is a normahzation factor. In case the correlator fulfills S{t) oc S{t), the process is 
a Gaussian process with a white noise. Rather, ri{t) is supposed to be a stationary process 
and the autocorrelation function S{t) refiects a colored noise behavior. That fact offers some 
profound differences with the conventional Langevin equation with a white noise |23] . As for 
instance the autocorrelation time tc of ri{t) is not zero, the stochastic process governed the 
evolution equation for x{t) is a non-Markovian one. Due to the correlator S{t) the values 
of i]{ti) for ti > to with i]{tQ) is fixed, depends on all the values of r]{t) inside the interval 
to <t < ti. Although our result is independent on the special realization of S(t), we consider 
to be specific the following colored noise behavior 

Sit) = -^expi-\t\/Q. (6) 

The system is assumed to be coupled to a heat bath with the fixed temperature T. 
For the subsequent discussion we need the Greens function related to Eq. (pi): 

G{t-t') = Q{t-t')g{t-t'), 
with g(t — t') = sinuJeU — t'), k = , Ue = \/uJn — k^ . (7) 

In terms of the Greens function the solution of Eq. ^ reads 

/oo 
dt'G{t-t')[7]{t') + K{t')]. (8) 

-oo 

The work performed on the system up to time t is given by 

W{t;T,Q= [ dt'K{t')x{t'). (9) 

^— oo 

Apparently the work depends on both the time r, compare Eq. (pi), and on the correlation 
time tc of the colored noise according to Eq. (|6| . 
Due to the initial condition ^(t = 0) = it results 

x{t) = f_ dt'j^g{t-t')[r^{t') + K{t')]. (10) 



We are interested in the work performed during the entire interval — oo < t < oo. Inserting 
Eq.(lO) into Eq. ^ the work function can be expressed by W{oo; r, tc) = Wri^r; t^) + Wk{'t) 
with 



Wk{t) 



d 



dt'K{t',r)J dt"-g{t'~t"Ut") 
'K{t',r)l' dt"^g{t'-t")Kit"] 



dt 



(11) 



Whereas the work due to the external force, denoted as Wk{t), depends only on the duration 
time r, see Eq. (Is]), the stochastic part Wrj{T,tc) is influenced by both time scales r and tc- 



In particular, the first relation of Eq. (11) offers a multiplicative coupling between the work 



protocol K{t; r) and the noise //(t). While in case of a white noise both quantities K{t) and 
ri{t) are instantaneously coupled at the time t, in the case of a colored noise the work is 
coupled to the whole history of the noise. Eventually this coupling leads to the breakdown of 



the JR for the system driven by the colored noise. In performing the integration in Eq. (11) 



we have to take into account the different time behavior of the external force according to 
Eq. (g. It results 



W^{t,Q 



R{t,r) 



T 



dtR{t; T)7]{t) with 

0(r-t) f dt'G{t' ~t). 

Juiax{o,t) 



(12) 



For the subsequent calculation the Fourier transformed function R{uj) is required which is 
given by 



R{u;) 



^o(e^"" - 1) 



K — l[U — Up 



K^ + iuj 



Up, 



K — i{uj + Ulf.) 



(13) 



Using the working protocol we can calculate the work function Wk{t) according to Eq. (11) 
The result is 



WAr) 



, ,4^2 



I -4 
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cos U^pT 
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(14) 



Now we are able to study in detail the behavior of the system. 



III. GENERALIZED JARZYNSKI RELATION 



According to Eq. ^ the total work is defined by 

/oo 
dtR{t)7]{t) + WK{r) . 
-oo 



Here the function -R(t), see Eq. (12), can be rewritten as 



/oo 
K{t')dt>G{t' - t)dt' . 

Following Eq. ([I]) let us consider the expression 

/g-/3H/(oo;r)^ ^ l^^-P Sr^dtR{t)-n{t)x ^-l3WK{r) 



(15) 



(16) 



(17) 



Using Eq. ([s]) the average in the last equation can be easily performed. To that aim we 
remark that after making Fourier transformation the exponential term in the probability 
measure Eq. ^ reads 



duo 

2^ 



Tj * {ijj)S{uj)ri{uj) . 



Performing the integration over the probability distribution we obtain 

(exp(-/5iy (oo; r, Q) = exp (^ j^ dcu ^^|?4^) exp(-/3iyx(r)) . 



(18) 



The correlation function for the colored noise, defined in Eq. ([6]), reads after Fourier trans- 
formation 



S{u;) 



h 



(19) 



1 + {ujt.Y ■ 

In the limiting case tc = the Gaussian noise behavior is retained. Using the last expression 



combined with Eq. (|13|) we get the contribution from the noise in the following form 

2 



duo I R{uo) 



27r S{uo) 



-^ {i + 2kt + (uootcY - 4 - 



- e"''^ cos{uOeT) [1 + {uoQtcY - 4 



UOo 



Ue 



3 + {uootcY - 4 
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LdQ 



(20) 



To discuss the JR let us introduce the function V{T,tc) according to 

(exp(-/5W^(oo;r,g)^exp(/?V^(r,g) with V(r, tj = /? T ^^^ 

Joo ^TT b{uj) 



WK(r). 
(21) 



In case the JR is fulfilled this quantity should be identified with the free energy V{t, tc 



-AF. Inserting Eqs. (20) and (14) into Eq. (21) it results 
V{T,tc) = ~-^^ 2( -1 \l-e- 



~2^ 



T 



K 






sin(cc;er) ] — 1 



(22) 



As usually assumed the system is in equilibrium before the external force K{t) is applied. 
After fixing the force K{t) = Kq for t > r the system is again in an equilibrium state. 
Therefore we can calculate the difference of the free energy. The partition function is given 
by 

Z{K) = ^ I dxdpe-^^ . (23) 

From here we get the difference of the free energy between the equilibrium initial and final 

state 



AF = (3-^ In 



ZiO) 



Kim 



(24) 



Notice that the sign of AF is chosen in such a manner that it corresponds to the free energy 
of the damped oscillator applied to the system from outside. To illustrate the behavior of 
our system let us introduce a new function U by 



f/(r,g = \/(r,g + AF. 



(25) 



If this function IJ is nonzero for any finite r and tc 7^ the JR has to be modified. Inserting 



Eq. (22) and Eq. (24) it results 



[/(r,tc 
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(26) 



The function f/(r, t^) > is a measure for the dissipative energy. In case the correlation 
time tc of the colored noise tends to zero we find U{t, t, = 0) = 0. That case corresponds 
to a white noise behavior and consequently the JR is retained. Likewise for an infinitely 
slow process r ^ cxo it results U{t — > 00, tc) = and the the JR is fulfilled again. In the 
last case the system remains in equilibrium. The behavior of the dissipation U is shown in 
Fig. ni Notice that the ratio ^ plays an decisive role for the breakdown of the JR. The JR is 
fulfilled for — —* 0, i.e. the characteristic memory time of the colored noise process should 
decay within the time of the work protocol. However in all other case, especially for a finite 
correlation time tc, the conventional JR relation is violated. As argued before in that case 
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FIG. 1: Function U according to Eq. (25) as function of kt for different ratios ^^ = 20 (solid); 10 
(dashed) and 5 (dotted line) . 

the work performed on the system becomes history dependent and a more general relation 
is fulfilled: 

For a non-Markovian process the conventional JR is modified by the correlation time of the 
colored noise. Let us stress again that our result is not restricted to the special choice of 
the colored noise due to Eq. ^ with a correlation time tc- From Eqs. (18) and (19) one 
concludes that for a Gaussian process S^u) is independent on u and consequently the JR is 
fulfilled. In case >S'(ti;) depends on the frequency uj the JR has to be modified accordingly. 



IV. CONCLUSIONS 

In this paper we have demonstrated that the conventional Jarzynski relation, which 
is broadly discussed in a variety of different fields, has to be corrected if the system is 
subjected to a colored noise. For such a realization the history dependent stochastic process 
offers a significant influence on the work performed on the system. Additionally to the fact 
that the work is not a state variable as already in the conventional statistical mechanics, 
the work function becomes even history dependent leading to additional contributions. 
Originated by the coupling of work function and noise the correction term to the JR is 
determined by the autocorrelation time of the colored noise. In view of our result the proof 
of the validity of the Jarzynski in |22] is restricted to a special class of systems with memory 
and not valid for all kinds of non-Markovian processes. Concerning that problem we are 
aware that there exits different interpretations. In the stronger classification proposed by 
van Kampen [23] systems with memory and white noise should be considered as Markovian. 
Here we have considered a simple model of non-interacting damped oscillators subjected to 
a colored noise which is doubtless a non-Markovian process. 

This work has been supported by the DFG (SFB 418). M.S. is grateful to the University 
Halle for kind hospitality. 
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